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ABSTRACT 

We present Chandra observations of the galaxy cluster AC 114. The cluster shows a strongly irregular 
morphology, with signs of multiple merging activity. We report the discovery of a soft X-ray filament 
originating close to the core of the cluster, curving to the south east for approximately 1.5'. We 
confirm that X-ray emission is associated with two of three mass concentrations identified in previous 
gravitational lensing studies of this object ijNataraian et aljll998tlCamnusano et aLlfeOOl^ . These two 
mass concentrations are located at opposite ends of the soft filament, and we interpret this as evidence 
for interaction between them. The northern part of the cluster reveals signs of further dynamical 
activity. Two sharp discontinuities are detected in both the surface brightness and temperature 
profiles, evincing another, more recent merger event which took place in, or close, to the cluster core. 
A preliminary combined mass and lensing analysis shows, in spite of the merger activity, remarkably 
good agreement between lensing and X-ray masses computed with the assumption of hydrostatic 
equilibrium. 

Subject headings: gravitational lensing - Galaxies: clusters: general - Galaxies: clusters: individual: 
AC 114 - X-rays: galaxies: clusters 



1. INTRODUCTION 

AC 11 4 (Abell S1077) i s a well-studied galaxy 
cluster l|Couch fe Sharpies! Il987|) : at a redshift 
z = 0.313, it has a re latively large velocity-dis persion 
a = 1660 km s" 1 (jMahdavi fc Gelled l200l and a 
compact core dominated by a cD galaxy. Its strong 
lensing power was first noticed by ISmail et"aT1 l|1991fl 
in a survey for bright gravitational lensing arcs. Sub- 
sequently, several br ight arcs and multiply i maged 

sources llSmail et all 119951 iNataraian et alJ Il998t 
iCamnusano et al.ll200lD were discovered. 
In this paper we analyze the first Chandra observation 
of AC 114. This hot (kT = 8.0 keV) cluster shows an 
irregular morphology with a highly elliptical shape. The 
central cD galaxy, slightly displaced with respect to the 
centroid of the X-ray emission, is aligned in the general 
direction of the X-ray brightness elongation. 
Throughout this paper we quote errors at the 90% 
confidence level and, unless otherwise stated, we use 
H Q = 72 km s^Mpc -1 (O m = 0.3, fl A = 0.7). This 
implies a linear scale of 4.46 kpc per arcsec at the cluster 
distance z = 0.313. 



2. OBSERVATION AND DATA ANALYSIS 

The galaxy cluster AC 114 was observed on 2000 
November 17 by the Chandra X-ray Observatory, us- 
ing the Advanced CCD Imaging Spectrometer (ACIS) 
(Observation-ID 1562) for 74.3 ksec; the cluster emission 
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lies mainly in the back-illuminated CCD S3. 
The Chandra Interactive Analysis of Observation (CIAO 
2.3) software was used throughout. Level 2 event files 
supplied by the Chandra X-ray Center (CXC) were up- 
dated with new gain maps; bad pixels were removed us- 
ing the observation-specific bad pixel list provided by the 
(CXC). Periods of high background rate were eliminated 
by visual inspection of the light curve in the energy range 
10.0 — 12.0 keV; time intervals with rates higher than 
70 cts/ (100 s) were eliminated, leading to a final useful 
exposure time of 71.7 ksec. 

Throughout the spectral analysis the blank field 
ba ckground sets ( aci5*5 D2000-01-29bkgrndN000*.fits) 
bv lMarkevitch et alJ l)2003ft available in the Chandra Cal- 
ibration Database (CALDB) were used, after projecting 
them onto our observation. The blank fields were taken 
at different locations and in different periods, so the fol- 
lowing two corrections were made: 

• To check for temporal changes in the particle back- 
ground rate, we compared the count-rates in the 
cluster and blank fields in the 10.0— 12.0keV energy 
band, where negligible cluster emission is expected. 
We found a count-rate ratio of 1.015, which, though 
remarkably close to unity, was nevertheless used to 
scale the blank field data. 

• We then looked for a possible residual background 
due to the soft (E < 2.0 keV) Galactic component 
of the cosmic X-ray diffuse emission which varies 
with position in the sky. We extracted spectra from 
the AC 114 field and from the blank fields inside a 
region where the lowest emission from the cluster 
was expected and no other cosmic source was ob- 
served. In this region (at ~ 0.8 Mpc from the clus- 
ter center) a count-rate c/r Q 3 7 = 0.179 cts s _1 
was measured in this region, very similar to the 
one measured in the blank fields (c/r 3 7 

I). 1 71 cts s _1 ), confirming that the X-ray emission 
in this area is mainly due to background. We found 
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Fig. 1. — Chandra image of AC 114 in the 0.3 - 7.0 keV energy 
band, binned to 2" pixels. 



an excess flux (~ 10%) in the cluster field, rela- 
tive to the background field, at energies less than 
2 keV. The difference between the two spectra was 
appropriately scaled and then subtracted from each 
source spectrum throughout the spectral analysis. 

3. CLUSTER MORPHOLOGY 

The Chandra image of AC 114 in the 0.3 - 7.0 keV en- 
ergy band is shown in Fig. ^ The emission arises mainly 
from the intra-cluster component and from a few addi- 
tional point sources. While the X-ray surface bright- 
ness decreases rapidly out of the cluster center toward 
the north, and especially the north-east direction (where 
at just a few arcseconds from the cluster center the X- 
ray emission decreases sharply almost to the background 
level), diffuse, elongated fainter emission stretches to- 
ward south/south-east down to the edge of the field of 
view. Since the cluster clearly shows a strongly irregular 
morphology, we first looked for any energy dependence. 
We created adaptively smoothed, background and expo- 
sure corrected images of the area covered by ACIS detec- 
tor S3 in three energy bands: E < 0.3 keV, 0.3 — 2.0 keV 
and 2.0 — 10.0 keV as shown in Fig. |3 (from top to bot- 
tom, respectively). Pixel values of all detected point 
sources were replaced with values interpolated from the 
surrounding background regions; the CIAO tools wavde- 
tect and dmfilth were used for this purpose. The images 
were then adaptively smoothed with minimum and max- 
imum signal-to- noise thresholds set at 3 a and 5 a. 

In the softest energy band [E < 0.3 keV) the emission 
is dominated by two extended components: the cluster 
and an elongated tail-shaped feature extending southeast 
from the cluster. The elongated southeastern component 
becomes less prominent at higher energies. 
The centroid of the overall X-ray emission (computed 
within a radius of 4' in the 0.3 — 10.0 keV band) is located 
at a = 22 h 58 m 48! 1, S = -34° 47' 59"4 (J2000.0), but 
the cluster does not show a single X-ray peak; notice- 
able emission is associated with the cluster cD galaxy, 
which is located about 10" south-east (22 h 58 m 48!4, 




Fig. 2. — Adaptively smoothed, exposure-corrected images of the 
central 5 arcmin of the cluster in three energy bands: E < 0.3 keV, 
0.3 — 2.0 keV and 2.0 — 10.0 keV (from top to bottom, respectively). 
In the top panel Clump 1 and Clump 2 show the positions of 
two additional clumps of matter as estimated with strong and weak 
lensing measurements (see § 171. 
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Fig. 3.— HST image of the cluster with 0.3 - 2.0 keV Chandra 
X-ray contours superposed in green. The positions of the cluster 
centroid and of the central cD galaxy are plotted as a blue circle 
and as a black box, respectively. White and blue crosses are at the 
locations of the X-ray peaks in the E < 0.3 keV and 0.3 — 10.0 keV 
energy bands, respectively. 



-34° 48' 08'.'5 J2000.0). The positions of the cluster 
centroid and of the central cD galaxy are plotted as a 
blue circle and as a black box, respectively, on the HST 
image of the cluster in Fig. [3] Superimposed in green are 
the contours of the X-ray emission in the energy band 
0.3-2.0 keV. The offset between the cD galaxy and the 
centroid of the diffuse intra-cluster medium (ICM) emis- 
sion is evidence that the latter is most probably not in a 
state of hydrostatic equilibrium. 

3.1. Two-dimensional model fitting 

To parameterize the peculiar morphology of AC 114, 
we modeled the emission in the central 4' of the cluster 
in the 0.3 — 10.0 keV band. All point sources in the field 
were removed. Using the SHERPA software we fitted 
the cluster surface brightness to several models. Sin- 
gle and d ouble, circular and elliptical, two- dimensional 
(3- models l|Cavaliere fc Fusco-Femianolll976j) have been 
used. These take the general form: 



S(x,y) = S 1 



sjx\ + yj/e 2 



-3/3+1/2 



with: 



Xi = (x — x c ) cos 9 + (y — y c ) sin9 
Vi = (y ~ Vc) cos 9 - (x - x c ) sin 9 

where r c is the core radius, (3 the slope parameter, e the 
axial ratio, 9 the position angle (north over east) and x c 
and y c the location of the center of the model. The X-ray 
emission turns out to be best described by a simultaneous 
fit of two separate components. The northern component 
is centered on the centroid of the cluster X-ray emission 
and has an axial ratio e = 0.61. The center of the south- 
ern component is located ~ 0.9' south-east of the X-ray 
centroid and is centered on the elongated soft tail visible 
south of the cluster. This southern component is more 



Fig. 4. — Left panel: 0.3 — 10.0 keV smoothed image with over- 
plotted the linearly spaced isointensity contours of the residuals af- 
ter subtracting the best-fit model of the northern component from 
the data. Right panel: Image of the best-fitted two-dimensional, 
double elliptical /3-model of the surface brightness; superimposed 
are the linearly spaced isointensity contours of the two components 
of the /3-model. 



elliptical (axial ratio e = 0.50) than the northern one, 
and is rotated about 15° with respect to it. The po- 
sitions, shapes, orientations and relative amplitudes of 
these two model components are listed in Tabled The 
sum of the components is shown in a false color image 
in the right panel in Fig. 01 superposed are the isointen- 
sity contours of each component. To better illustrate the 
shape of the softer and fainter southern component, we 
subtracted the best-fit model of the northern component 
from the data in the 0.3— 10.0 keV band. Smoothed con- 
tours of the residuals are shown in the left panel in Fig. 01 
over-plotted on the smoothed data. The regular, ellipti- 
cal shape of the southern component is evident. For the 
remainder of this paper this Southern component will be 
called the "Tail" . 



4. SPECTRAL ANALYSIS 

4.1. Global Cluster Properties 

A global spectral analysis for the cluster was per- 
formed inside a circle centered at a = 22 h 58 m 48? 1, 
S = -34° 47' 59'.'8 (J2000.0) and within a radius of 
2.5'; this region includes both the cluster and the 
tail. The emission from all detected point sources 
inside this region was excluded. We fitted the data 
to an a bsorbed XSPEC is otherm al plasma emission 
code by iKaastra k, Mewd 11993T) including the FeL 
calculations of lLie<iaTii*et al.l (| 1995) 5 , folded through 
the appropriate response matrices, corrected for the 
ACIS time dependent absorption due to molecular 
contamination. An excess of soft emission is observed in 
the spectrum at energies below 0.5 keV, most probably 
due to the emission of the soft Tail. We were able to 
obtain a good spectral fit only when energies below 
0.5 keV were excluded. (We return to this point in 
Section 14.31 below.) Results obtained in this spectral 
analysis are summarized in Table [21 For the spectral 
range E > 0.5 ke V, results a re in agreement with 
previous findings by lAlTenl (|2000^ . 

The hydrogen column density along the line of sight is 
not well constrained by our fit, and may be subject to 

5 McKaL code hereafter. 
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Table 1. Fit parameters of the 2D elliptical double- f> model 



Component 




X, 


Vc 




r c 





e 


e 


A 








(J2000.0) 




['] 






n 


[cts s~ 1 


C)- 2 ] 


Northern 


a = 


= 22 h : 


58 m : 


48? 1 


1.0 ± 0.8 


1.7 ± 1.6 


0.61 ± 0.08 


167 ± 7 


0.19 ± 


0.03 




5 = 


= -34° 


: 47' 


: 59'.' 8 














Southern 


a - 


= 22 h : 


58 m : 


50? 3 


1.3 ± 1.0 


0.6 ± 0.4 


0.50 ± 0.04 


143 ± 3 


0.086 ± 


0.008 




6 = 


: -34° 


: 48' 


: 46? 1 















Note. — Fit parameters of the 2D elliptical doublc-/3 model of the cluster: x ci y c is the central position of each component, r c the core radius, 
{3 the slope parameter, e the axial ratio, 9 the orientation angle (north over east) and A the amplitude of the model at x c ,y c . 



systematic errors in the calibration of the (changing) 
low-energy response of ACIS. We therefore analyzed 
the cluster spectrum in the restricted energy range 
1.5 — 7.0 keV, where the value of hh does not affect the 
fit. Results, shown in Tabled are consistent, within the 
uncertainties, with the fits over the 0.5 — 7.0 keV band. 
We conclude that spectral fits over the broader band are 
reliable. 



4.2. Flux and Luminosity 

The cluster emission can be traced out to a radius of 
3.5' ~ 940 kpc. Within this radius the cluster bolo- 
metric flux and luminosity (at the measured tempera- 
ture of kT = 8.0 ± 0.5 keV) are: /x(bol) = 5.6loj x 
10~ 12 erg s- 1 cm- 2 , L x (bol) = 12.9±^ x 10 44 erg s" 1 . 
AC 114 follows the typical Lx-T cluster rela- 



tion, with an expected L^ ip (bol) = 12.7 ± 2.6 x 
10 44 erg s _1 l)Arnaud fc EvrarcTl 999'l in agreement with 
our findings. In contrast, we find that AC 114 does 
not fall on established trends relating X-ray proper- 
ties to velocity dispersion. If we adopt the higher 
vel ocity dispersion, a = 16 601}^ km s _1 , computed 
by iMahdavi fe Gellerl ((2001^ within a radius of 1 Mpc, 
bolometric luminosities and temperatures as high as 
L^P(bol) = 45.61 2 .- 4 x 10 4 4 ergs" 1 and kT exp = 
13 ± 5 keV a re predicted <|Mahdavi fc Gellerl 12001 
IWu et all Il999|) . Even the lower velocity dispersion 
(a = 1380 1? 28 km s" 1 ; compu t ed wi thin 0.9 Mpc) re- 
ported by IGirardi fc MezzettH l)2001t) leads to an ex- 
pected value of Lx P ( bo1 ) = x 1q44 er S s_1 > sig- 
nificantly higher than the value we measure; the expected 
gas temperature of kT exp = 10 ± 4 keV is instead consis- 
tent to the measured value. 

The mean Lx — <? and a — T relations presumably 
characterize dynamically relaxed clusters, so departures 
from these relations m ay be expected for unrelaxed sys- 
tems fW u et aLlll999yi . Our results are thus consistent 
with the view that AC 114 is not fully relaxed. While 
the significant discrepancy between the published veloc- 
ity dispersions might suggest that this conclusion ten- 
tative at best, may also, in itself, be indicative of the 
dynamical complexity of the cluster. 

4.3. The Tail 

In order to characterize the spectral properties of the 
elongated Tail, we computed and then compared spec- 
tra extracted within two elliptical regions: one centered 
on the Tail and the second located at 22 h 58 m 44!0, 



—34° 47' lO'.'l, diametrically opposite the Tail with re- 
spect to the cluster center. Both regions have shapes and 
orientations given by the Southern component in Tabled 
the semi-major axis is equal to 0.5' for each region. The 
two resulting spectra, with the respective best-fit MeKaL 
models superimposed, are shown in Fig.0 the spectrum 
extracted from the Tail is plotted as a black line, while 
the one from the region in the north is shown in red. 
The MeKaL model was fitted in the 0.5 — 7.0 keV energy 
band. Resulting best-fit parameters for both regions are 
consistent with the global cluster spectrum (see Table|2J. 
Although the two regions are symmetric respect to the 
cluster centroid, their emission spectra differ significantly 
(see Table |21 and Fig. |3J). First, the Tail is brighter by a 
factor of three. Second, the Tail shows a strong excess 
of emission with respect to the fitted model at energies 
below 0.4 keV; no such soft excess is present in the north- 
ern region. 

In the Tail region, we have a superposition of the emis- 
sion from the Tail and from the cluster. At energies above 
1.5 keV, the spectra of the two regions are similar, but 
the Tail is brighter. This result suggests the Tail region 
contains plasma with a temperature similar to that of 
the Northern region. It is thus natural to conclude that 
the ICM in these two regions are dynamically related, 
and that the cluster emission is more elongated in the 
southern direction, possibly due to some recent dynamic 
interaction. 

In order to characterize the excess soft emission in the 
southern region, we subtracted the two spectra extracted 
above, after having rescaled the northern component to 
have the same flux as the southern component at ener- 
gies higher than 1.5 keV. The limited number of photon 
counts, and calibration uncertainties in the very low en- 
ergy range, make a precise characterization of the spec- 
tral properties of the Tail rather difficult. An estimate of 
its temperature, however, has been obtained from broad- 
band colors. Color ratios among three energy bands 
(0.1 - 0.8 keV, 0.8 - 1.5 keV and 1.5 - 8.0 keV) were 
computed for simulated spectra as a function of temper- 
ature. We find a direct correspondence between the color 
ratio CR = counts[(0.8-1.5) keV]/counts[(l. 5-8.0) keV] 
and the temperature 6 . The observed CR = 2.97 implies 
a temperature for the Tail of kT = 0.7 keV and a bolo- 
metric luminosity, within 0.5' of 2.7 x 10 43 erg s _1 . We 
stress that the given values of the temperature and of 
the bolometric luminosity of the Tail are only tentative 
estimates since an exact spectral characterization of the 

6 kT = EXP (1.003 ± 0.003 - 2.48 ± 0.04 x log(CR)). 
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Table 2. MeKaL spectral fits for the cluster central 2.5' and for two ellitpi- 
cal regions, one centered on the tail (region a) and the second one symmetric 
respect to the cluster centroid (region b). 



Region 


Energy 


kT 


Metallicity 


Rcdshift 


HI column density x 2 /d.o.f. 




[keV] 


[keV] 


[Z ] 




[I0 20 cm- 2 ] 



0.313 (fixed) 0± 2 292/271 

0.34+°;°2 0± 2 282/270 

0.313 (fixed) 1.31 (fixed) 343/272 

0.313 (fixed) 7 ± 10 177/203 

0.313 (fixed) 4.1 ±0.6 187/93 

0.313 (fixed) 1.31 (fixed) 221/94 
0.313 (fixed) 3.9 ±1.2 31/39 

0.313 (fixed) 1.31 (fixed) 39/40 



CI. center 0.5-7.0 8.0 ± 0.5 0.25 ± 0.09 

CI. center 0.5-7.0 8.0 ± 0.5 0.30 ± 0.09 

CI. center 0.5-7.0 ll.OtJg 0.20 ±0.12 

CI. center 1.5-7.0 8.61, 1 ,- 2 0.24 ±0.10 

A 0.1-7.0 7.5lg | 0.37 ±0.20 

A 0.1-7.0 12.212 2 0.37 ±0.33 

B 0.1-7.0 13±5° 0At° l 

B 0.1 - 7.0 26l^ OA^ll 



region could not be obtained. Furthermore normaliz- 
ing the spectra extracted in the two different regions and 
creating count-rate ratios using two energy ranges can be 
affected by systematic errors depending on uncertainties 
in the background estimate which can affect the North- 
ern and the Southern regions in different ways and on 
the choice of the threshold temperatures. 
A soft filament of similar shape, but of smaller dimen- 
sions (~ 8 kpc) was observ ed in the galaxy cluster 
Abell 1795 ijFabian et al.ll2001[) and explained by the au- 
thors as a possible cooling wake, caused by the motion 
of the central cD galaxy. In the case of AC 114 the high 
bolometric luminosity of the Tail, and the fact that it ex- 
tends for ~ 400 kpc, tends to exclude the motion of the 
cD galaxy as a possible cause of its origin. An inf ailing 
group, stripped of its intra-group medium when it passed 
close to the cluster core is a more likely explanation for 
the soft component in this cluster. 

Spatially dependent degradation in the quantum effi- 
ciency and calibration uncertainties could in principle be 
affecting our imaging and spectroscopic results at such 
low energies. Our detection of a soft excess is spatially 
confined in the region of the Tail; significant sensitiv- 
ity fluctuations at that specific position of the instru- 
ment would be required to cause such an effect, which 
should therefore be detectable also in other observations. 
We have analyzed spatial photon distribution in a low 
(E < 0.3 keV) and in a broad energy range (0.3 — 7.0 keV) 
both in the blank fields and in a deep observations of 
a galaxy cluster. In both cases we did not detect flux 
increase at the location of the Tail. Spectra extracted 
from the blank fields within the two elliptical regions 
(described at the beginning of this section) also show no 
significant difference at low energies. Calibrations issues 
and spatially dependent degradation in the instrument 
quantum efficiency should therefore not be affecting our 
detection of the soft Tail. 



5. A RECENT MERGER CLOSE TO THE CLUSTER CORE? 

Fig. n shows a sharp cutoff in the cluster brightness 
toward the north/north-east, while the cluster emission 
is more diffuse and extended in the other directions sug- 
gesting that the cluster, or an internal substructure, is 
moving toward the north/north-east direction creating a 
sharp edge in the front and a more diffuse feature behind. 
To explore this possibility, we extracted surface bright- 



Energy [keV] 

Fig. 5. — Spectra, with superimposed best-fit models (top panel), 
and residuals (bottom panel), extracted from the region in the 
south where the Tail is located (black line), and from a region in 
the north symmetric respect to the cluster center (red line). 



Table 3. Count-rates in different energy bands 



Region 


c/r (0.1 


- 0.4 keV) 


c/r (0.4 


- 8.0 keV) 




[cts 


sec -1 ] 


[cts 


scc~ 1 ] 


south 


30.9 


x 10~ 4 


4.0 


< 10~ 2 


north 


7.5 


x 10~ 4 


1.5 


< 10~ 2 



ness profiles from three contiguous wedges in the north- 
ern portion of the cluster. The wedges are plotted on 
panel e) of Fig. superposed on the 0.3 — 2.0 keV 
smoothed image of the cluster; the annuli are centered 
at the cluster X-ray centroid and have an axial ratio 
b/a = 0.8. The resulting surface brightness profiles, plot- 
ted in red, light blue and black in panel a) of Fig. 
are extracted from the north-west, north and north-east 
wedges, respectively. Although all three profiles show 
a steep decrease as a function of radius, the north-east 
profile reveals two discontinuities: a sudden decrease of 
a factor of rs 4 at a distance of around 0.1 Mpc from the 
cluster center, and a second smaller decrease between 0.3 
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and 0.4 Mpc. 

In the area between these two jumps the temperature of 
the ICM shows a maximum (panel c) in Fig.|HJl. Although 
the errors in the temperature profile are large, the tem- 
perature increase can be clearly seen as the sharp hard 
curved front in the hardness ratio (HR) map (panel f) 
of Fig. EJ); the two dashed black lines show the positions 
where the two discontinuities, both in surface brightness 
and temperature, are observed. A second hard region 
is observed on the western side of the cluster core, but 
neither a sharp edge in the HR map, nor drastic jumps 
in the surface brightness and temperature, are observed 
there. The HR map has been obtained from the ratio 
of the adaptively smoothed images in the 1.0 — 10.0 keV 
band and below 0.5 keV, with point sources excluded; the 
use of a smoothing technique with a fixed sigma does not 
change features in the final HR map. In the resulting im- 
age, colors from dark blue to red and yellow, correspond 
to increasingly harder ratios. 

Similar discontinuities in the surface brightness, with 
corresponding jumps in the ga s temperature, have 
been observed in othe r clusters ijVikhlinin et alJ |2001t 
iMarkevitch et all 120021) . Also, recent adiab atic hydro- 
dynamical simulations bv lMathis et alJ l)2003|) show how, 
during a cluster merger, two shock fronts form, and sub- 
sequently expand from the central regions into the ICM, 
right after the cores of two merging structures pass one 
another; AC 114 shows very similar features, with two 
hard regions located on opposite sides of the X-ray cen- 
troid. While the western hard region is not as well de- 
fined, in the eastern one we have seen how a sharp front 
with a corresponding temperature increase is clearly vis- 
ible. The two discontinuities in the surface brightness 
emission of the cluster are observed in projection on the 
plane of the sky. To obtain an estimate of the posi- 
tions and amplitudes of the two jumps, we approximate 
the density distribution in the cluster with an analytical 
model. We use a model composed of three power laws 
centered in the cluster centroid, having different ampli- 
tudes and slopes in three regions: interior to the inner 
discontinuity, between the two discontinuities, and exte- 
rior to the outer discontinuity. We then fit the X-ray 
surface brightness distribution predicted by this density 
model to the observed profile. The free parameters are 
the three power law indexes, the two radii and the two 
amplitudes of the jumps (a±, a 2 , oes, Ti 2 , r23, Ai and 
A2 respectively). The best- fit surface-brightness model, 
which provides a reasonably good fit to the data, and 
the associated (3D) density model, are shown in pan- 
els a) and b) of Fig. H3 The best-fit radial slopes are: 
at = 0.90 ± 0.19, a 2 = 0.89 ± 0.25 and a 3 = 1.80 ± 0.15 
and the positions of the two jumps are at distances 
i'i2 = 0.108 Mpc and r 23 = 0.401 Mpc from the cluster 
center. The amplitude of the jumps in the gas density 
are easily computed from the square root of the relative 
amplitudes of the three best-fit surface brightness mod- 
els: 

— = 2.1 ±0.5; — = 1.49 ±0.16 

92 P3 

At ri2 we observe a sharp gas density decrease and an 
increase in the gas temperature which, together with the 
only small decrease observed in the pressure, indicates 
that ri2 does not correspond to a shock front. 
At r 2 3 we instead observe a sharp decrease of both the 



gas density and the temperature. This is consistent with 
the physical conditions expected if the gas was experi- 
encing a shock. 

AC 114 could therefore be going through a merging pro- 
cess in its core, with consequences sim ilar to what is ob- 
serve d in clusters such as Abell 2142 l)Markevitch et all 
2000). Within r i2 we would then be observing the core 
of one of the two sub-clusters, moving toward the north- 
east, after a relatively recent merger. If this sub-cluster is 
moving supersonically, a bow shock would form in front 
of it, with a subsequent sudden decrease both in surface 
brightness and temperature, as we observe at ^3. 
Given the relatively large errors in the temperature and 
in the fit of the projected gas density model, it is im- 
possible to draw precise quantitative conclusions. We 
have, nevertheless, made crude estimates of the shock pa- 
rameters. From the density jump across the bow shock 
at r23 we estimate the Mach number of the core to be 
M = 1.3 ± 0.1. This should correspond to a temperature 
and pressure ratio across r23 of AfcTff* = 1.3 ± 0.2 keV 
and APJf* = 1-8 ±0.4 erg cm -3 . The observed tempera- 
ture and pressure discontinuities across r23 are AkT^ s — 
1.4±° j keV and AP° 3 bs = 2.0±[Jjerg cm" 3 which are con- 
sistent with the estimated Mach number. 
A second method f o r com puting the Mach num- 
be r due to iMoeckel lll949fi and applied to clusters 
bv lVikhlinin et alJ ^OOlfT vields a similar estimate: M = 
1.15. 

6. MASS DISTRIBUTION 

The irregular morphology and the apparently unre- 
laxed dynamical state of AC114 suggest that the ICM 
may not be in hydrostatic equilibrium. We therefore es- 
timate the mass profile in three different ways. We begin 
with the simple assumption that all the observed hot gas 
is part of a single structure (Global model), first mod- 
eling its surface brightness distribution with a single [3- 
model and subsequently with a double /3-model. Finally, 
we consider the hypothesis that the cluster is composed 
of two dynamically independent systems -the Northern 
component and the Southern Tail- (N+S model). 
For the first step of our mass analysis, we extracted ID 
surface brightness profiles within elliptical concentric an- 
nuli spaced by 5 arc-seconds, with ellipticity, position 
angle and center matching those of the northern com- 
ponent (see Table QJ. The distance r from the clus- 
ter center is measured along the major axis of the el- 
lipses. The profiles were extracted in the 0.3-7.0 keV en- 
ergy range. We extracted background profiles from the 
blank field data, in the same regions and energy band, 
scaled as described in Section above. Corresponding 
profiles extracted from the exposure maps provided the 
required average values of the effective exposure time 
within each annulus. The best-fit single /3-model to 
the resulting surface brightness profile is described by: 
S = 0.253 ±0.016 [cts s -1 arcmin -2 ], /? = 0.423 ±0.007, 
r c = 0.37 ± 0.03 arcmin, x r 2 cd . = 4 - 3 - A double (3- 
model increases the quality of the fit, providing an accu- 
rate description of the cluster profile despite its strong 
irregularity: S01 — 0.213 ± 0.016 [cts s -1 arcmin -2 ], 
01 = 0.9 ± 0.5, r cl = 0.557 ± 0.045 arcmin, S 02 = 
0.056 ± 0.008 [cts s -1 arcmin -2 ], (3 2 = 0.60 ± 0.06, 
r c2 = 1.8 ± 0.4 arcmin, x^cd = 1-1. 
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Fig. 6. — a): surface brightness profiles extracted in the north-west, north and north-east wedges (plotted in red, light blue and black, 
respectively; superimposed as a solid black line is the best-fit projected model of the north-east profile), b): density model, c): temperature 
profile extracted in the north-west wedge, d): temperature profile extracted in a 150° wide wedge (centered at the cluster X-ray centroid) 
in the southern portion of the cluster, e): adaptively smoothed color image of the cluster in the energy band 0.3-2.0 keV; superposed are 
the wedges within which the three surface brightness profiles, shown in the top left panel, have been extracted. The pink curves show the 
position of the fronts ri2 and T23- f): HR image where colors from blue, to red an yellow indicate increasingly harder ratios, as labeled in 
the color scale at the bottom of the image. The two fronts rj2 and T23 are plotted as black dashed black curves. 
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The electron density n c is then obtained by inverting 
the model which describes the surface brightness distri- 
bution. In the spherically symmetric case t he inversion 
of the double /3 model has been resolved bv lXue k, Wul 
(2000). Since we are dealing with elliptical distributions, 
we had to extend the formalism to a more general ellip- 
soidal symmetry; we considered an oblate ellipsoid with 
the z-axis along the line of sight, and we did not con- 
sider a possible inclination of the principal axes. Under 
these assumptions, the electronic density for the single 
/3- model reads 



1 a i \ I m 

n e p (m) = n e0 1 + 



where m = 



% + z 1 defines an ellipsoid coordi- 



nate. For a double concentric /3-model 



1/2 



,20 



(m) = < ^2 [nl f3 (m;n e oi,r ci ,Pi)]' 



At large radii from the cluster center (~ rsoo: radius 
corresponding to an overdensity A = 500 with respect 
to the critical density) an assumption of a prolate ge- 
ometry leads to negligible differences in the cluster total 
mass and gas mass fraction; differences in the total mass 
of clusters, even in compressed or elong ated geometries, 
have been estimated to be less than 4% (jPiffaretti et al.l 
2003) . Once the gas distribution is known, we can obtain 
the total mass density, under the hypothesis of hydro- 
static equilibrium. Towards North-East the ICM shows 
a sudden temperature increase between 0.2 and 0.4 Mpc 
probably due to a recent merger (see § EJ- Such phe- 
nomenon does not appear to affect the Southern portion 
of the cluster, which shows instead a much more regu- 
lar temperature profile (see panel d) in Fig.^J; we have 
therefore approximated the overall ICM as isothermal. 
The total dynamical mass density corresponding to each 
independent, isothermal gas phase is obtained from: 



Ptot 



\ 4nGfim p ) 



V 2 In) 



where ks is the Boltzmann constant, G the gravitational 
constant, T gas the intra-cluster gas temperature and fxm p 
the mean particle mass of the gas. For the single and 
double /3-model, using the cluster average temperature 
(Table [21, we obtained (within a sphere with 1 Mpc ra- 
dius) the following values of the total and gas masses: 

Global, 1/3 : M tot [l Mpc] = (3.7 ± 0.3) • lO 14 ^ 1 M Q 
M gas [l Mpc] = (7.4 ±0.5) • I0 13 h^ /2 M© 

Global, 2/3 : M to t[l Mpc] = (4.0 ± 1.7) • lO 14 /^ 1 M© 

M gas [l Mpc] = (7.9 ± 1.6) • lO 13 ^ 572 M© 

If the assumption of isothcrmality is dropped and the 
average cluster temperature profile is used, the resulting 
cluster total mass is consistent within 10% with the above 
values for radii larger than 250 kpc. At smaller distances 
to the cluster centre the observed temperature is lower 
than the cluster average temperature, leading to a sys- 
tematically lower cluster total mass (within 25% of the 
above values, for radii larger than 75 Mpc). The use of a 



temperature profile averaged over 360°, causes the tem- 
perature increase observed in the North-East, which is a 
spatially localized phenomenon, to be instead associated 
to the whole cluster, leading nevertheless to systematic 
errors. The assumption of isothermality is therefore used 
for the remainder of this paper. 

Finally, we introduced the Tail as a separate component 
in our mass analysis. We modeled the cluster as two dy- 
namically independent systems, with different electron 
temperatures, seen in projection along the line of sight. 
These are strict assumptions, but only within this frame 
can we investigate the effect of the elongated Tail on 
the mass distribution of the system. Surface brightness 
profiles were extracted within two sets of annuli shaped 
according to the results of 2D-fit (see Table [TJ. The pro- 
file for the Northern component was extracted within a 
180° centered toward the North direction, where the con- 
tribute of the Tail is negligible. A corrected ID profile of 
the Tail was obtained subtracting the profile of the north- 
ern component, taken into account their relative shift and 
orientation. Best-fit parameters of the resulting surface 
brightness profile of the Northern component are for a 
single-/? model: Sqn — 0.32 ± 0.03 [cts s _1 arcmin -2 ], 
N = 0.389 ± 0.007, r cN = 0.226 ± 0.026 arcmin, 
X 2 cd . = 4.0 and for a double-/? model: S ni = 0.24 ± 
0.02 [cts s" 1 arcmin -2 ], f3 N1 = 1.3 ±0.3, r cN1 = 0.84 ± 
0.19 arcmin, S 0N2 = 0.029 ± 0.004 [cts s -1 arcmin -2 ], 
j3 N2 = 0.65 ± 0.18, r cN2 = 3.1 ± 1.1 arcmin, x 2 cd . = L6 - 
The Tail is instead best described by a single-/? model: 
S os = 0.068 ± 0.005 [cts s" 1 arcmin -2 ], S = 1.8 ± 1.0, 
r c s = 2.8 ± 1.0 arcmin, Xi- C d = 0.8. Temperatures of 
8.0 keV and 0.7 keV were used for the Northern and 
the Southern component, respectively. Under these as- 
sumptions, the total mass density was computed adding 
linearly the mass of the tail to that of the Northern com- 
ponent (each one obtained independently), taking into 
account their relative positions and orientations. While 
in the inner regions around the X-ray centroid the mass 
density is dominated by the Northern component, at 
large radii the contribution from the tail becomes sig- 
nificant. This is due to both the relative positions of the 
two components and to the smooth profile of the Tail in 
contrast to the sharply decreasing one of the Northern 
component. The resulting values of the total and gas 
masses are: 

N + S, 1 /? : M tot [l Mpc] = (4.5 ± 1.1) • 10 14 h^ 1 M© 



M gas [l Mpc] 



.4 ±2.6) ■ 10 Li h 



131,-5/2 



72 



Mr; 



Consistent values of the masses are obtained using a 
double-/? model for the Northern component. Henceforth 
with N±S we will imply N+S, 1 model. 
Although these mass estimates are highly model- 
dependent, some general qualitative trends emerge. 
When we assume that all the emission comes from a sin- 
gle structure (Global model), the results of the mass ob- 
tained with a single or a double /?-model are similar, es- 
pecially at small radii (Fig. 0) ; this reflects the fact that 
in the central region of the cluster a simple 1-/? model 
still provides an acceptable approximation to the surface 
brightness profile. 

If we assume instead that the system is a superposition 
of two distinct systems, we obtain systematically higher 
values of the mass, especially in the inner regions. This 
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quite general result is a consequence of the non-linear de- 
pendence of X-ray luminosity on cluster masses. As seen 
for the Global model, also for the Northern component 
both single and double /3-models provide a good descrip- 
tion of the surface brightness profile leading to similar 
final values of the masses. 

We compute the gas mass fraction at r 50 o; to ease 
comparison with previous work, we quote results for 
H = 50 km s^Mpc -1 (fi m = 0.3, Q A = 0.7). We 
remark that going out to rsoo needs a radial extrapola- 
tion of ~ 30%. For the three models described above 
wc find: ^ haX ' 10 (r 5OO ) = 0.33 ± 0.04, fgf^rsoo) = 
0.31 ± 0.09, fg+ s (r 50 o) = 0.39 ± 0.14. The values of 
r5oo for the three models are 1363 kpc, 1499 kpc and 
1474 kpc, respectively. Considerably lower values of 
fp-asfrsn n) ~ 0-15 ± 0.05 are estimated for relaxed clus- 
ters bv lEttori et"aD l|2003ft . 

7. COMPARISON OF X-RAY AND STRONG- AND 
WEAK-LENSING MASS ESTIMATES 

From the results of a spectroscop ic survey of faint 
lensed galaxies in the core of AC 114. iCampusano et alJ 
(2001) revised the lensing mo del of the cluster de- 
rived bv lNataraian et all l)1998[) which uses both strong 
and weak lensing co nstraints. The revised mass re- 
construction made bv ICampusano et al.l l|200lD includes 
the strong constraints determined from the redshifts 
of two multiple-image systems: the 3-image S (z = 
1.8 6), that had been alr eady included in the model 
bv lNataraian et all lll998j). and the newly discov ered 5- 
image E (z = 3.347). ICampusano et all l)2001|) found 
that a mass distribution composed of an elliptical cluster- 
scale component (the central cluster) plus an additional 
bimodal component, with the two clumps centered on 
two other galaxy concentrations, reproduces the observed 
lensed images. The locations of these two hypothetical 
additional clumps of matter are shown in the top panel in 
Fig. |21 While we detect no excess X-ray emission above 
the background from Clump 2, diffuse X-ray emission 
can definitively be seen at Clump 1 . It is particularly in- 
teresting to note that Clump 1 is located at the Southern 
edge of the elongated soft Tail. 

7.1. Comparing Mass Estimates 

From their detailed analysis of weak and s trong lens- 
ing data for AC 1 14, iNataraian et al.l l)1998j) estimated 
the cluster mass within 75.0, 150.0 and 500.0 kpc of the 
cluster center. Their results are listed in Table0|together 
with our values for the projected X-ray masses. 
The X-ray and lensing masses are in good agreement, 
for the two-component X-ray models, even at the small- 
est radii. In particular, ratios Mi cns /Mx at 75 kpc are 
1.23 ± 0.16, 1.1 ± 0.5, and 0.95 ± 0.16 for the Global-1/3, 
Global-2/3, and N+S models, respectively. The results are 
remarkable in view of results from some other unrelaxed 
clusters in which the X-ray-determined masses are only 
one-thir d to one-half the lensing masses at radii l ess than 
100 kpc HMa.chac.ek et al.ll2002jlXue fc Wull200l. 
The use of a detailed map of the surface mass density is 
cruc ial to accurate modeling of strong lensing phenom- 
ena. ICampusano et ail (|200 If ) used the galaxy density to 
predict a cluster mass distribution that reproduced the 
observed strong-lensing images. 
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Fig. 7. — Mass profiles computed assuming the Global-l/3-model, 
Global-2/3-model and the N+S model, plotted as a dotted, dashed 
and solid line, respectively. The 90% confidence levels for the 
N+S model are plotted as d ot-dashed lines. Results from lensing 
from Nataraian ct al. (1998) arc also plotted. 



Table 4. Lensing and X-ray mass measurements 



Radius 


M x 


M x 


M x 


M !c „ s 


[kpc] 


Global 1-/3 


Global 2-/3 


N+S 




75.0 


3.4 ± 0.3 


3.7 ± 1.5 


4.4 ± 0.6 


4.2 ± 0.4 


150.0 


9.7 ± 0.7 


10 + 4 


11.2 ± 1.7 


12 ± 1.5 


500.0 


39 ± 3 


41 ± 16 


47 ± 10 


40 + 4 



Note. — Column 1: radius. Cols. 2-4: Projected X-ray masses. 
Column 5: Mi» nfl from lNataraian et al.lll99gD . All results in this Table 
are given for: Hq — 50 km s — Mpc - . All values of the masses are 
given in [10 13 Mq]. 



Since the X-ray emission supplies an independent tracer 
of the underlying cluster mass, we used our inferred (X- 
ray-derived) mass distribution for a new strong lensing 
analysi s. For this p urpose we used the gravlens soft- 
ware l)Keetonll200l() . We found that the mass distri- 
bution computed from the Global model is sub-critical, 
and is therefore unable to produce any strong lens- 
ing. The N+S model produces two lips-shaped caustics, 
one in side the other, orient ed perpendicular to one an- 
other l)Schneider et al.l ll9921 as shown in Fig. IH1 Sources 
inside the inner caustic are multiple imaged in patterns 
similar to the ones observed in the E-systems. A better 
match with the lensed images is obtained when the cen- 
ter of the Northern component is slightly shifted toward 
the direction of the cD galaxy; this provides a reason- 
ably good reproduction of both E and S systems with 
a resulting average distance between the lensed images 
and their best-fit positions: <Ie = 2.4" for the five images 
in the E-system and ds = 2.9" for the three images in 
the S-system (see Fig. |SJ . A detailed new strong lens- 
ing analysis of AC 114, based on the X-ray mass results 
obta ined in this paper is the obj ect of a forthcoming pa- 
per ijSereno fc De Filippisll200^ . 

' http://astro.uchicago.edu/ckeeton/gravlens 
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Fig. 8. — Lensing properties of the N+S model, for the two 
sources at z = 3.347 (5-images E system) and z = 1.86 (3-images 
S system). The black and gray lines represent critical curves and 
caustics, respectively. Crosses show the image positions of the E- 
system, X the ones of the S-system; filled stars and triangles the 
corresponding image positions predicted by the N+S model. The x- 
axis coincides with the North. Distances are in units of arcseconds. 
The origin coincides with the position of the X-ray centroid. 



8. SUMMARY AND CONCLUSIONS 

We have presented the first Chandra observation of 
the galaxy cluster AC 114, which has revealed a strongly 
irregular X-ray morphology. Below 0.5 keV the X- 
ray emission is dominated by two main components: 
the cluster, roughly centered on the optical position of 
AC 114, and a diffuse " Tail", extending almost 400 kpc, 
from the cluster center to the south-east. The Tail has 
an energy distribution different from the rest of the clus- 
ter, showing significant excess of soft emission. The Tail 
connects the cluster core with the location of Clump 1, 
an additional clump of matter previously inferred from 
lensing analysis. We also detect diffuse X-ray emission 
from Clump 1, confirming its existence. We conclude 
that Clump 1 is probably associated with the cluster, 
and may have interacted with it. 

We propose a possible scenario in which Clump 1, lo- 
cated nearby a large over-density of galaxies, in its mo- 
tion from the north-west through the cluster, has been 
ram-pressure stripped of most of its intra-group gas, now 
still visible as the soft Southern Tail. During its interac- 
tion with the cluster, Clump 1 would have also distorted 
the intra-cluster gas, causing the asymmetrical stretch of 
the cluster emission toward south-east. Such a picture is 
support ed by recent adap tive mesh hydro/N-body simu- 
lations ijBurns et al.ll2003|) which show soft extended fea- 
tures similar to our tail forming after substructure merg- 
ers. 

AC 114 also shows further signs of recent dynamical ac- 
tivity close to the cluster center. A sharp decrease in the 



cluster surface brightness toward north-east and a corre- 
sponding peak in the gas temperature are observed, sug- 
gesting the motion of a sub-structure toward this direc- 
tion which compresses and heats intra-cluster gas ahead 
of its path. 

These two putative interactions take place along differ- 
ent directions in the plane of the sky. Their time-scales 
are also different. High-temperature fronts forming in 
front of fast moving sub-structures are short-living phe- 
nomena, and are therefore signs of "rec ent" merging, 
proces ses; comparison with simulations bv lMathis et alJ 
(2003) lead to an estimated time of the merger of ~ 
200 Myr. In the Tail region there is instead no sign of 
high velocities involved, and considering the length of the 
Tail, the closest approach of Clump 1 to the core must 
have been a much earlier event (« 800 Myr ago if one 
assumes the velocity of Clump 1 to be « 500 km s~ ). 
These two merger events were thus likely to have been 
independent phenomena. We note that in a hierarchical 
cosmological scenario, in which clusters form from ac- 
cretion of matter through sets of intersecting filaments, 
multiple merging processes, directed along different axes, 
are expected. 

In spite of this evidence for complex interactions, we 
have obtained estimates of the X-ray masses that are 
in remarkably good agreement with lensing results, even 
at distances as small as 75 kpc from the cluster center. 
The simplest isothermal model for the X-ray emission is 
under-critical and thus cannot reproduce the lensing re- 
sults. Two-component models, in which the components 
interact, or are seen in projection along the line of sight, 
however, provide a consistent estimate of the total mass 
and approximately reproduce the positions of strongly 
lensed background sources, though the best agreement 
is found if the X-ray emitting gas is displaced slightly 
with respect to the potential minimum. Agreement be- 
tween X-ray and lensing mass estimates has recently been 
fo und also in another i rregu lar cluster (MS1008. 1-1224) 
bv lEttori Ez Lombardil l|2003|) . We therefore advocate for 
the robustness of the X-ray mass estimates, and conclude 
that the assumption of hydrostatic equilibrium can yield 
accurate mass estimates even in clusters as dynamically 
active as AC 114 once the gas density distribution is 
properly mapped. 
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